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Abstract. Nucleus-nucleus interaction potentials in heavy-ion fusion reactions are extracted from 
the microscopic time-dependent Hartree-Fock theory. When the center-of-mass energy is much 
higher than the Coulomb barrier energy, extracted potentials identify with the frozen density approx- 
imation. As the center-of-mass energy decreases to the Coulomb barrier energy, potentials become 
energy dependent. This dependence indicates dynamical reorganization of internal degrees of free- 
dom and leads to a reduction of the "apparent" barrier. Including this effect leads to the Coulomb 
barrier energy very close to experimental one. Aspects of one-body energy dissipation extracted 
from the mean-field theory are discussed. 
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INTRODUCTION 

The interplay between nuclear structure and dynamical effects is crucial to properly 
describing fusion reactions at energies close to the Coulomb barrier. Coupled-channels 
models [1,2] have been widely used to describe the entrance channel of fusion reactions. 
While in general rather successful, these models have in common several drawbacks. 
First, nuclear structure and dynamical effects should be treated in a unified framework. 
Second, important channels should be guessed a priori. Mean-field theories based on 
the Skyrme energy density functional provide a rather unique tool for describing nuclear 
structure and nuclear reactions in a unified framework, i.e., all of the dynamical coupling 
effects between collective and intrinsic degrees of freedom. In nuclear reactions, appli- 
cation of the time-dependent Hartree-Fock (TDHF) to heavy-ion fusion reactions was a 
major step |j3iH,|5|,|6|]. Since recent TDHF calculations can now include all terms of the 
Skyrme energy density functional used in static Hartree-Fock calculations n7l.[8l.l9l.[l0l]. 
the description of nuclear reactions using TDHF should be revisited. 

In this contribution, as one of the applications of TDHF to nuclear reactions, we 
propose a method for simultaneously extracting nucleus-nucleus potentials and friction 
coefficients associated with one-body energy dissipation from the microscopic TDHF 
theory [LL]. In this method, we assume that fusion dynamics is described by one- 
dimensional macroscopic dissipative dynamics on the relative distance between collid- 
ing nuclei. In order to validate our assumption, we first compare extracted potential 
with alternative mean-field methods |T2, Hfl. Then, we discuss the property of nucleus- 



nucleus potential and one-body energy dissipation deduced from TDHF. 



METHOD 



The potential and friction coefficient of one-body energy dissipation are extracted as 
follows: (i) The TDHF equation of head-on collision is solved to obtain the time evo- 
lution of the total density of colliding nuclei, (ii) After dividing the total density into 
two densities at the separation plane defined in Ref. [H], we compute at each time dif- 
ferent macroscopic variables: relative distance R, associated momentum P, and reduced 
mass /i. (iii) We assume that the time evolutions of R and P obey a classical equation of 
motion including a friction term which depends on the velocity R: 

dR P dP dV 

where V{R) and y{R) denote the nucleus-nucleus potential and friction coefficient, 
respectively. The friction coefficient y{R) describes the effect of energy dissipation from 
the macroscopic degrees of freedom to the microscopic ones, (iv) Equation ([T]) has 
two unknown quantities dV/dR and y{R). These quantities are obtained by using two 
TDHF evolutions with slightly different energies llll]. The potential V{R) is deduced 
by integration over 7? using its asymptotic Coulomb potential at large relative distances. 
For the TDHF calculations, the three-dimensional TDHF code developed by P. Bonche 
and coworkers with the SLy4d Skyrme effective force [7] is used. The mesh sizes in 
space and in time are 0.8 fm and 0.45 fm/c, respectively. This method is called hereafter 
dissipative-dynamics TDHF (DD-TDHF). 



NUCLEUS-NUCLEUS POTENTIAL 
Comparison of potential 

The present method assumes that the mean-field dynamics can properly be reduced 
to a one-dimensional macroscopic dissipative equation. In order to validate this assump- 
tion, we first compare the potential that we obtained, denoted by V^^, with other tech- 
niques based on mean-field theories. 

The potential V^^{R) is displayed by the solid line in Fig[T]for the ^^0-1-^^0 reaction 
at the center-of-mass energy ^c.m. = 34 MeV. The potentials obtained by the density- 
constrained TDHF (DC-TDHF) method lll2ll (dashed line) and by the frozen-density 
(FD) approximation [13] (filled circles-dotted line) are also shown for comparison. In 
DC-TDHF, dynamical effects are partially accounted for by minimizing at each time step 
the total energy under the constraint of the density reached along the TDHF path. The 
FD approximation is based on the sudden approximation and estimates potential energy 
from energy density functional with the condition that projectile and target densities 
are frozen to their respective ground state densities at each R. Figure [T] shows that the 
potentials extracted from DD-TDHF and from DC-TDHF are almost identical even well 
inside the Coulomb barrier. This gives confidence in the specific macroscopic equation 
[Eq. ([U)] retained to reduce the microscopic dynamics. In addition, both methods are 
almost identical to the FD approximation (for R > 6.5 fm). This indicates that little 
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FIGURE 1. Comparison of potential energies for the '^O+'^O reaction obtained from different models. 
The solid, dashed, and filled circles-dotted lines co rresp ond to our result, to the density-constrained 
method lfl2ll . and to the frozen density approximation lfl3ll . respectively. 



reorganization of densities occurs in the approaching phase at this energy (^c.m. = 
34 MeV), which is well above the Coulomb barrier energy. As a consequence, the 
Coulomb barrier predicted by TDHF is almost identical to the one obtained in the FD 
case (the difference being less than 0.1 MeV). It is worth mentioning that our method 
assumes neither sudden nor adiabatic approximation. 



Energy dependence of extracted potential 

To illustrate the center-of-mass energy dependence of the potential, Fig. |2] presents 
potentials obtained with DD-TDHF using several center-of-mass energies ranging from 
Ec.m. = 55 MeV to 100 MeV for the "^'^Ca-^'^'^Ca reaction. Again, in the high energy 
limit, potentials identify with the FD case. In addition, an increase of center-of-mass 
energy from ^c.m. = 90 to 100 MeV leads to identical results indicating the stability of 
DD-TDHF as the energy increases. In opposite, as ^c.m. decreases, potentials deduced 
from DD-TDHF deviates from the FD case. As ^c.m. approaches the Coulomb barrier 
energy, a small change in ^c.m. significantly affects extracted potential as illustrated by 
the two energies ^c.m. = 55 MeV and 57 MeV shown in Fig. [21 

This effect is a direct consequence of reorganization of densities in the approach- 
ing phase. This is clearly illustrated in Fig. [3] where density profiles obtained for the 
^*^Ca+^^Ca reaction at ^c.m. = 55 and 90 MeV are shown for specific R values. In Fig.[3l 
only the case of ^c.m. = 90 MeV resembles the FD case, remaining spherical. At low en- 
ergy ^c.m. = 55 MeV, a clear deviation from the FD profile is observed. As the two 
partners approach, deformation of the two nuclei takes place. This deformation initiates 
the formation of a neck at larger relative distances compared to ^c.m. = 90 MeV. This 
center-of-mass energy dependence of the extracted potential reflects the difference in the 
density profiles accessed dynamically during the mean-field evolution. Note that similar 
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FIGURE 2. Potential energy for the ^"Ca+^^Ca reaction extracted at different center-of-mass energies. 
The FD potential is shown by the filled circles-dotted line. 
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FIGURE 3. Density profiles obtained from TDHF for different relative distances R = 10.26 (left), 9.82 
(middle), and 9.52 fm (right) for the '^^Ca+'^^'Ca reaction at ^c.m. = 55 (top) and 90 MeV (bottom). 



dependence is a priori also expected in the DC-TDHF method [[l2l[l4ll . which accounts 
for the dynamical deformation of the densities. 

Dynamical effect on extracted potentials is systematically found in all reactions con- 
sidered here. Figure |4] shows the difference between the barrier height deduced from 
DD-TDHF and the barrier height from experiment [15] as a function of extracted bar- 
rier height for the i^O+^^Ca, '^^^'^^Csi+'^^^'^^Csi, and ^^Ca+^OZr reactions. 
The solid line is the result when potential barrier is extracted in the high energy limit 
of DD-TDHF (^c.m. ^ Vg), whereas the dashed line for the low energy limit of DD- 
TDHF (£c. m. ~ Vb)- Dynamical reduction of the barrier energy is clearly seen for all the 
reactions. Moreover, the value of the barrier energy approaches the experimental data. 
This underlines the importance of dynamical effects close to the Coulomb barrier and 
illustrates the degree of precision of our technique. 
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FIGURE 4. Difference between the barrier height deduced from DD-TDHF and experimental barrier 
height ifTsll as a function of extracted barrier height for the reactions indicated in the figure. Vb is deduced 
from high energy TDHF (solid line) and from low energy TDHF (dashed line). 



FRICTION COEFFICIENT FROM MICROSCOPIC MEAN-FIELD 

In most practical models with dissipation applied to nuclear reactions at energies around 
the Coulomb barrier, the mechanism of energy dissipation is assumed to be of one-body 
type, where energy dissipation is caused by collision of nucleons with the wall of mean- 
field potential and by nucleon exchange between colliding nuclei, i.e., the so-called wall- 
and- window formula fl^ [itI \lS]. TDHF includes the mechanism of one-body energy 
dissipation from the microscopic point of view because of the self-consistency of mean- 
field. Therefore, we investigate the property of energy dissipation from the microscopic 
point of view by DD-TDHF [19]. 

In Fig. m we present reduced friction parameters defined as P{R) = y{R)/ii{R) as 
a function of R scaled by the Coulomb barrier radius Rg. As the colliding nuclei ap- 
proach, the magnitude of the friction coefficients monotonically increases. Figure [5] 
clearly shows that the order of magnitude of /3(i?) and the radial dependence are al- 
most independent on the size and asymmetry of the system. Besides, we compare our 
results with a microscopic model [ 20 ] based on the linear response theory by the fiUed- 
circles. They agree very well. DD-TDHF gives energy dissipation in reasonable order of 
magnitude and points out that extracted friction coefficients have universal behavior. 



SUMMARY 

A novel method (DD-TDHF) based on the macroscopic reduction of TDHF has been 
used to extract nucleus-nucleus potential as well as one-body energy dissipation. The 
DD-TDHF gives important insight in the dynamical effects and provides a way to extract 
friction coefficient from dynamical microscopic theory. 
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FIGURE 5. Extracted reduced friction parameter I5{R) — yiR) I P^{R) as a function of R scaled by the 
Coulomb barrier radius Rg for several reactions. A microscopic friction by Adamian et al. |20] is shown 
by the filled-circles for comparison. 
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